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SPECIAL  FINITE  ELEMENTS  FOR 
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by 

K.  C.  WATTERS 


''V  SUMMARY 

This  paper  describes  a  special  hybrid  finite  element  In  represent  a  leaded  hale  in  a 
sheet.  The  shape  functions  of  the  element  satisfy  stress  equilibrium  and  strain  compati¬ 
bility  throughout  the  element,  and  the  applied  loading  boundary  conditions,  l  hc  applied 
loading  is  represented  by  a  finite  Fourier  series,  anti  the  coeffit  ictus  of  the  clement  shape 
functions  are  matched  with  the  Fourier  coefficients.  )  computer  program  lias  been  written 
to  generate  the  element  stiffness  and  stress  recovery  matrices,  llte  program  also  produces 
equivalent  nodal  forces  and  initial  stresses  related  to  the  applied  loading.  Act  urate  results 
have  been  obtained  from  several  example  analyses. 
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1.  INTRODUCTION 


In  this  paper,  a  special  finite  element  enclosing  a  loaded  hole  in  a  plane  sheet  is  described  in 
some  detail.  This  work  had  its  origins  in  an  attempt  to  develop  a  special  element  to  enclose  a 
fragment  damage  site  in  an  aircraft  structure:  such  an  clement  would  facilitate  residual  strength 
analyses  of  damaged  aircraft  structures  for  use  in  vulnerability  assessments.  For  example,  vs  hen 
an  aircraft  is  attacked  by  certain  types  of  warhead  the  damage  incurred  by  the  aircraft  skin 
consists  of  many  closely  spaced  fragment  perforations  w  liich  can  be  idealised  as  holes  with  cracks 
radiating  from  them.  A  natural  first  step  was  to  consider  the  case  where  the  damage  consists  of 
uncracked  circular  holes  in  a  plane  sheet.  Since  problems  involving  sheets  with  circular  holes, 
both  unloaded  and  loaded,  are  of  common  occurrence  in  aircraft  structures,  quite  apart  from  the 
particular  application  which  gave  rise  to  the  present  study,  these  were  considered  in  de'ail:  see 
also  Reference  I. 

The  basic  concept  of  the  finite  clement  method  of  structural  analysis  is  that  a  complex 
structure  is  partitioned  into  a  network  of  simple  shaped,  standard  elements.  The  stillness 
characteristic  of  the  whole  structure  is  then  formed  by  assembling  the  known  stiffness  charac¬ 
teristics  of  the  individual  elements;  see  Reference  2.  However,  for  regions  of  high  stress 
gradients,  such  as  occur  around  holes,  it  has  proved  beneficial  to  develop  special  element  types, 
rather  than  to  use  fine  meshes  of  standard  elements.  Rao  era/.'’ 1  were  among  the  earliest  to  apply 
this  concept.  In  Reference  3.  a  general  method  for  formulating  special  elements  around  stress 
concentrators,  including  stress  singularities,  is  presented.  For  this  method,  the  displacement 
state  within  the  element  is  defined  by  functions  which  satisfy  conditions  of  stress  equilibrium, 
strain  compatibility,  and  any  boundary  conditions  (BCx)  of  the  stress  concentrator.  Such 
elements  arc  termed  “hybrid  elements"  in  Reference  3.  in  that  they  combine  concepts  of  con¬ 
tinuum  mechanics  with  finite  elements.  A  consequence  of  the  use  of  rather  complex  displacement 
functions  for  hybrid  elements  is  that  boundary  displacements  are  incompatible  with  those  of 
adjoining  standard  elements.  However,  a  comparison1  of  two  sector  elements,  one  with"natural 
mode"  shape  (i.e.  displacement)  functions  and  the  other  with  simple  conforming  shape  functions, 
has  clearly  demonstrated  the  superior  performance  of  the  former  despite  its  violation  of  inter¬ 
element  displacement  compatibility. 

In  Reference  3.  Rao  el  al .  demonstrate  the  use  of  a  special  hybrid  element  surrounding  an 
unloaded  circular  hole  in  plane  stress  problems.  In  the  present  paper  that  work  is  extended  to 
the  case  of  loaded  circular  holes.  The  paper  is  set  out  as  follows.  A  general  outline  of  hybrid 
clement  theory  is  given  in  Section  2.  Then,  in  Section  3.  the  detailed  application  of  that  theory 
to  derive  an  element  for  a  loaded  rir<-ular  hole  is  given,  and  in  Section  4.  the  associated  computer 
program  is  full  described.  Resu'i  _i  sample  runs  of  the  program  are  presented  in  Section  5.  and 
a  general  discussion  is  given  in  Section  6. 


2.  HYBRID  ELEMENT  THEORY 

For  the  stress  analysis  of  two-dimensional  elastic  structures,  the  governing  equation  of 
equilibrium  and  compatibility  is:’ 

V'^=0  (2.1) 

where  <f>  is  the  stress  function  and  V  is  a  differential  operator,  defined  in  polar  co-ordinates  by: 

I  >  I  !>  - 

r  .V  /■-  ,su- 


1 


(2.2) 


\  separated  variable  solution  is  proposed  for  <jt  in  the  form: 


where 


i  =  F(r)G{«) 


(2.3) 


0(0)  = 


|  sin  nO | 
jeos  nOy 


(2.4) 


The  upper  and  lower  functions  in  the  bracketed  notation  are  each  separately  valid.  On  solving 
(2.1).  (2.3)  becomes: 


in  =  [o»r"  ;-6nr  »  +  c„r-  »  -  dnr-  (2-5) 

with  the  two  degenerate  cases: 

(fin  =  do  r/>olnr  :  fo'-’e  </»r-lnr  (2.6) 

ii  =  [air  i  bir1  t  ftr3  t  dtlnr)|^"  ^|.  (2.7) 

The  general  frequency  parameter,  n.  is  unrestricted  at  this  stage.  In  Reference  5  are  listed  several 
other  solutions  to  (2.1),  not  in  the  form  of  (2.4),  which  together  with  (2.5),  (2.6)  and  (2.7),  make 
up  the  full  general  solution  to  (2.1)  in  polar  co-ordinates.  These  other  solutions  are  grouped 
below  as: 


it  =  (’or'-0r/o<H  eirrtj  (2.8) 

Using  standard  formulae,5  the  stresses,  strains  and  displacements  can  be  derived  from  the  stress 
function.  Having  derived  these  quantities,  the  boundary  conditions  of  a  stress  concentrator  can 
be  satisfied.  For  example,  an  unloaded  hole  has  the  boundary  conditions; 

<?rr  =  r:u  =0  at  r  =  ro  all  0.  (2.9) 

These  conditions  allow  two  of  the  arbitrary  coefficients  a„,  b„.  <  „  and  ii„  to  be  determined  in 
terms  of  the  other  two.  The  allowable  values  for  n  are  restricted  to  integers  only,  to  satisfy  a 
requiremen  for  single-valucdness  in  a  multiply-connected  region.  Conditions  of  symmetry  can 
further  restrict  the  allowable  values  for  n.  Once  these  boundary,  single-valucdness  and  symmetry 
conditions  have  been  applied  to  the  general  solution  to  (2.1).  a  set  of  functions  remains  which 
satisfies  equilibrium,  compatibility,  and  all  conditions  of  the  stress  concentrator.  For  the  general 
case  of  this  set  of  functions,  tti>.  radial  and  tangential  displacements  take  the  respective  forms: 
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|cos  nOj 


(2.10) 


r  =  [fi\bn  ■  fad,,] 


j  -cos  Il0\ 

|  sin  oft)' 


The  equations  can  be  grouped  in  matrix  form  as: 
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FIG  1:  GENERALISED  SPECIAL  ELEMENT 


I  lie  displacement  formulation  (2.15)  applies  to  all  points  within  the  element,  including  the  nodal 
points.  Suppose  (2.15)  is  applied  to  each  of  the  ,V  nodes  in  turn  by  substituting  the  appropriate 
co-ordinates  into  the  matrix  L.  A  vector  of  nodal  point  displacements  can  then  be  formed  as: 


(say ). 


(2. IX) 


If  the  number  of  arbitrary  parameters  chosen  to  form  a  is  equal  to  the  number  of  nodal  point 
displacement  degrees  of  freedom,  then  the  C  matrix  will  be  square.  This  enables  the  relation¬ 
ship  (2.18)  to  be  inserted. 

a  =  C  '8.  (2.19) 


In  order  to  assemble  this  special  element  w  ith  the  other  elements  forming  the  structure,  it  is  neces¬ 
sary  to  derive  a  stiffness  relationship  between  the  nodal  forces  q  and  the  nodal  displacements  8. 

q  =  AS.  (2.20) 

The  stiffness  matrix  A  can  be  found  by  minimising  the  potential  energy  of  the  element,  using  the 
principle  of  virtual  work. 

A  =  (C  ')r  j  HTDHd\C  1  (2.21) 

where  the  integration  is  over  the  volume  of  the  element. 

The  relationship  (2.20)  can  be  readily  transformed  from  local  polar  co-ordinates  to  global  carte¬ 
sian  co-ordinates  by  using  a  transformation  matrix.  T. 

q  =  TQ:  8  =  7A.  (2.22) 

Q  =  7  rA  7 A  (2.25) 

^  AA  (say).  (2.24) 

Once  the  stiffness  matrix  A  has  been  generated  for  a  special  clement,  it  can  be  assembled  with  the 
other  elements  comprising  the  structure  and  a  solution  found  for  the  displacements  A  using  a 
finite  clement  program.  The  stresses  can  be  calculated  at  any  point  within  the  clement  as: 

a  =  Dlla  =  777/C  >7A  =  .SA  (say).  (2.25) 


3.  LOADED  HOLE  ELEMENT 

3.1  General 

A  special  element,  to  surround  a  circular  hole  with  applied  loading  on  its  boundary  ,  is  described 
here  (sec  also  Ref.  I).  The  element  is  able  to  account  accurately  for  applied  loading  by  directly 
matching  applied  load  and  boundary  stress  distributions.  The  applied  load  distribution  is 
represented  as  a  finite  Fourier  series.  Stress  function  coefficients  of  the  element  are  related  to  the 
Fourier  coefficients  of  the  applied  loading.  The  stress  function  formed  from  this  relationship 
contains  arbitrary  parameters  which  lead  to  the  formation  of  the  clement  stillness  matrix,  and 
fixed  constants  which  lead  to  the  formation  of  initial  stresses  and  equivalent  nodal  forces  to  the 
applied  loading.  As  would  be  expected,  the  stiffness  matrix  so  formed  is  identical  to  that  of  an 
equivalent  unloaded  hole  element.  The  influence  of  the  applied  loading  is  therefore  expressed 
through  the  equivalent  nodal  forces  and  the  initial  stresses.  The  equivalent  nodal  forces  are  input 
loads  for  the  finite  clement  analysis,  and  the  initial  stresses  are  added  to  those  calculated  by  the 
finite  element  analysis. 

3.2  Element  Shape 

The  basic  shape  of  the  element  is  an  annulus  with  a  polygonal  outer  boundary  whose  corners 
(nodes)  lie  on  a  circle  concentric  with  the  hole.  The  number  of  nodes  and  their  angular  co¬ 
ordinates  can  he  varied,  as  can  the  aspect  ratio  of  the  circumscribing  circle  radius  to  the  hole 
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radius.  Three  variants  of  this  general  shape  are  considered,  according  to  the  symmetry  of  the 
problem.  For  double  symmetry  only  one  quadrant  is  required,  for  single  symmetry  a  scmi- 
annulus.  and  for  the  unsymmetric  case  the  full  annulus.  The  three  element  variants  are  sketched 
in  Figure  2.  Also  in  Figure  2.  the  boundary  conditions  corresponding  to  the  states  of  sv  nimetry 
are  shown.  These  boundary  conditions  actually  alter  the  number  of  degrees  of  freedom  of  the 
element  and  that  etTect  is  also  shown.  The  rigid  body  motions  applicable  to  the  three  element 
variants  are  shown  in  Figure  2. 


3.3  Fourier  Analysis 

The  applied  loading  takes  the  form  of  normal  and  shear  stresses.  n.A(H)  and  ra(iq.  applied 
to  the  hole  surface,  r  =  r„.  The  applied  loads  are  represented  by  finite  Fourier  series  as: 

\s 

o.n(h)  =  <Tao/ 2  •  V  {(Jan  COS  //ft  ira„  $'\n  llH‘(.  (3.1) 

n  I 
Ss 

Tail))  =  (ran1  2)  •  Y.  I  r„„  sin  ll»  Ta„  COS  lltl  \ .  (3.2) 

it  *  I 

The  form  of  (3.2)  is  slightly  different  from  the  usual  form  of  a  finite  Fourier  series.  (3. 1 ).  to  simplify 
the  imposition  of  boundary  conditions  in  Section  3.4  below.  The  series  cutoff.  A'*,  is  selected 
large  enough  to  obtain  an  adequate  representation  of  nail1)  and  Ta(<')- 


3.4  Boundary  Conditions 

The  radial  and  shear  stresses,  derived3  from  the  general  stress  function  of  Equation 
(2.5),  are  listed  below. 

nrT  =  [  a„(n-  n)r"  -  h„in-  ■  n)r  "  -  <•„(/;-  n  2 )r" 


</«(//-  •  n 


2),  "j 


|  sin  ni'\ 
|cos  no 


(3.3) 


Tj,t  =  {«„(/!-  -  ;/)/•"  -  h„in-  ■  II)I-  "  -  ■  (  ,,(//-  •  11)1"  - 


(liiin- 


")'■  "] 


cos  nli\ 
sin  ni ij ' 


(3.4) 


The  single-valuedness  requirement  for  a  multiply-connected  region  restricts  n  to  integer  values. 
The  boundary  conditions  for  the  loaded  hole  are: 


<ijr  =  nalD)  and  t,„  =  ra((')  at  r  =  m.  all  0. 


(3.5) 


These  boundary  conditions  are  enforced  by  equating  Fourier  coefficients  of  the  applied  loading 
to  stress  function  coefficients. 


iiuin-  n)ru"  -  />„(/!-  •  n)rt,  "  -  c„in-  ii  2)r«" 

</»(«-  •  II  2)r„  "  =  (Jan. 
On  in-  n)rii"  -  />,,(»-  ■  i/)ru  "  -  ■  <■„(/!-  •  //)/(," 

linin’  ll)lt I  "  =  Tan. 


Solving  for  an  and  c„. 


On  = 


ill  2  )ra„  •  lla an 
2 ill-  n)rt,"  - 


in  ■  I  loi  nil, 


(3  M 


(3.7) 

(3.S) 


n  *  Ta n 

2 (n  ■  I  )r«" 


nrn 


-h„  ■  (ii  I 


(3.9) 
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i)  Doubly  symmetric  (DS) 


(b)  Singly  symmetric  (SS) 


(c)  Unsymmetric  (US) 


2N-2  degrees  of  freedom 
v  =  0  at  0  =  0,^  all  r 
No  rigid  body  motion  allowed 


y 


Local  cartesian  axes 


x 


2N-2  degrees  of  freedom 

v  =  0  at  0  =  0,ir  all  r 

x  translation  rigid  body 
motion  only 


2N  degrees  of  freedom 
v  unrestricted 
All  rigid  body  motions  - 
x  and  y  translations  and  rotation 


FIG  2:  LOADED  HOLE  ELEMENT  VARIANTS 
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Note  that  parameters  </»  and  <  „  contain  constant  terms  related  to  the  applied  loading  hut  their 
relationships  to  parameters  />„  and  i/„  are  independent  of  the  applied  loading.  Let  the  constant 
terms  of  parameters  a„  and  <•„  be  labelled  .r,„  and  sj„  respect  ixely  : 


(//  2)ra/j  •  //f7aH 

(3.10) 

2 (n-  n)i»"  - 

w  •  Ta„ 

2</i  ■  1 )/,," 

(3.1  1) 

In  Equations  (3.6)  to  (3.11).  the  boundary  conditions  (3.5)  hate  been  imposed  for  the  lower 
trigonometric  functions  of  Equations  (3.3)  and  (3.4).  Let  the  parameters  associated  with  the 
upper  trigonometric  functions  of  (3.3)  and  (3.4)  be  distinguished  by  a  dash  superscript  . ,/..  u„  , 
bn  .  i’ii  •  ‘In  ■  The  boundary  conditions  (3.5)  can  be  imposed  in  a  similar  manner  for  these  para¬ 
meters.  and  identical  equations  to  (3.6)  to  (3.11)  result,  except  that  all  parameters  base  dash 
superscripts,  i.e.  a„  .  b„  .  <■„  .  il„  .  n;u,  .  r;,„  .  \|„  .  .  The  relationships  (3.H)  and  (3.9)  can  be 

substituted  into  Equation  (2.5)  for  the  general  stress  function. 

'A»  =  I  [.sin'"  ■  .ve„/-  "]  [  (a  ■  Ib'n  •  r 

[  iii\i  -r"  (n  I  )/ii  "  •  r- 

i.c. 

'I’ll  =  if>u II  '  <AiiH. 

Where  ,A„i>  is  a  fixed  stress  function  related  to  the  applied  loading,  and  i/i„n  contains  the  arbitrary 
parameters  b„  and  d».  and  is  the  same  stress  function  as  for  an  unloaded  hole.  Stresses,  strains 
and  displacements  can  he  derixed  from  </>„. 


mu  -r-  "]/>„ 
|  sin  nil | 


Icos  n"J‘ 


(3.12) 


(3.13) 


Special  cases  for  „  =  0  and  I  can  be  derixed  in  the  same  manner  as  ahoxo.  using  the  stress 
functions  of  (2.6).  (2.7)  and  (2.8).  Stngle-xaluedness  requires: 


Then 


and 


Solvinti. 


=  r/o  =  0;  c\  -  tl\. 

1  r 

(3.14) 

f/»o  =  [</<»  •  />oIn/-  •  n»r-  •  fat1]. 

(3.15) 

<Jrr  -  [/>()/•  “  •  2oi]. 

(3.16) 

'm  —  [  fur  ']• 

(3.17) 

-  •  2(  0  -  <r;,o  2. 

(3.18) 

./( i/*n  "  =  Tail  2. 

(3.19) 

i  n  =  (t:1o  *4  0  •  5/  n  -ho. 

(3.20) 

./ n  =  O'  5/‘ojT;,o  , 

(3.21) 

V|n  —  0'5/'o“7-aii  ;  vjn  —  <T:lo  4. 

(3.22) 

'/>(»<>  =  [V|o^  •  s-»ur-]. 

(3.23) 

c/>„n  -  [In/-  0 ■  5(r  /■<> )-]/>(>. 

(3.24, 

,  ,  .,  ,  ,  ,f  'in  ")  .  2 in  )  cos ") 

b\r  1  •  i  ir’  ■  f/|t  lnr  <  ■  i/,  •  \ 

(cost'  |  r  sin") 

(3.25) 

3.5  Arbitral-)  Parameters 

li  remains  to  c hoose  values  oi  .  t,  tor  in  ilic  u-itor  a  of  the  .irhitr;ir>  parameters  />„  and  </„. 
The  rule  followed  is  to  choose  consecutive  allowable  rallies  of  n.  ranging  from  the  lowest,  until 
all  the  available  degrees  of  freedom  ot  the  element  have  been  filled.  The  allowable  values  of  n 
and  the  applicability  of  the  upper  and  lower  trigonometric  functions  (dashed  and  undashed 
parameters)  are  determined  In  conditions  ol  symmetry.  Symmetry  conditions  can  he  expressed 
as  additional  boundary  conditions  I ItC'C  are  listed  below  for  the  three  loaded  hole  element 
variants  (see  1  ig.  2). 


uiisv  ill  me  I  r  iv. 

no  restriction. 

(3.33) 

siriL’ls 

symmetric  / 

(1  at  "  0.  rr.  all  r. 

(3.34) 

ilouhk 

sMiimetrie  / 

0  at  ''  0.  .1  ?r.  all  r. 

(3.35) 

These  boundarv  conditions  atVect  the  n-.mher  of  degrees  of  freedom  of  the  element.  Normally, 
if  an  e  emenl  has  A  nodes  it  has  2  V  degrees  of  freedom,  corresponding  to  u  and  r  displacements 
at  each  node.  However.  (3.34)  and  (3.35)  suppress  the  r  displacement  at  the  first  and  last  nodes, 
leaving  the  singly  and  doubly  symmetric  variants  with  only  2\  2  degrees  of  freedom.  Conse¬ 

quently,  when  (he  ('matrix  is  evaluated  according  to  (2. IK),  the  r  displacement  is  not  included 
in  di  and  d.v.  Finally,  the  applicable  rigid  body  motion  components  have  to  be  included.  All  these 
factors  have  to  be  accounted  for  in  forming  a  and  the  result  is  shown  in  Table  I  for  elements 
with  A  nodes.  The  notation  tor  the  rigid  body  motion  components  is  shown  in  Table  2. 

3.6  Fixed  Parameters 

It  is  convenient  to  assemble  the  lived  paiameters.  ,vi„  and  .vj„.  related  to  the  applied  loading, 
into  a  vector  an  (say),  in  the  same  manner  as  the  arbitrary  parameters  are  assembled  into  a.  l  or 
the  general  case,  from  (3.12)  and  (3.13)  we  have: 


“  •  Vj "] 

|  sin  iih | 

(3.36) 

|vos  n('| 

M |  sin  | 

' 1 " r  |cos  i:1*) 

,  n|sin/tt'| 

Icosr/t'j 

(3.3?) 

-  A'l/i/l  nU'H.n)  V| 

li  1 1  «  -  V ,,/j „  -Vj„  tin 

(3  38) 

IV, „  j 

^  n  '  flu 

1 1 1 

(3.34) 

U,„  1 

l  v..„  ) 

-  /  (sa>). 

(3.401 
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TABLE 


Loaded  Hole  Element  Parameter. 


Doubly 

Singly 

Lnsymmciric 

symmetric 

symmetric 

Additional  boundary  con- 

i=0 

<■  =  o 

None 

ditions 

at 

i '  =  0. 

I  77 

ill 

4=0 

77 

all  r 

all  r 

Appropriate  functions  for 

CON  2 Hi* 

COS  //ft 

sin  (/,  •  1  )4 

6:  n  =  0.  1.  2 _ 

cos  //ft 

Applicable  rigid  body  motions 

None 

A  only 

All 

Element  degrees  of  freedom 

2  A  2 

2.V  2 

2.V 

r  , 

bn 

f  A  1 

r/>„" 

hi 

a 

A 

A  odd 

/->' 

A  even 

hi 

</.. 

’ 

/>j 

</, 

/■a 

</> 

</:« 

hi 

,  t 

i  *  l 

\ 

A 

i 

b<  a  i)  j 

h.v  m 

< 

> 

d,  .\  n  j 

i/.V  2  1 

*< i, 

r- 

h.v  2  1 

A. S  1 

f/(  \  1 1  2 

</a*  2  1 

A  \  1 

A 

h  A  2 

l*-  '  =d 

/’A  2 

A 

h.\  2 

,/s  2 

L  A  J 

A 

/>A  I 

A 

1 

L  J 

L  A  J 

TABLE  2 


Rigid  Bod\  .Motion  Parameters 


Rigid  body  motion 

Notation 

u 

/' 

,v  translation 

A 

A  cos  4 

A  sin  4 

r  translation 

A 

A  sin  4 

ft 2  COS 

Rotation 

0 

_ 

Ar 

(3.41) 

(3.42) 


Now.  if  the  total  stress  function  matching  (he  applied  loading  is  written  as: 

then 

*/>»  —  — /  ,.nO„o  -  /  iiOu. 


I 


I 


) 
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I'hsplacemcnts.  strains  and  stresses  can  he  derived  fiom  (5.42)  a-. 

d„  -  /.,, ,a„. 

c, i  :  //, ia„.  (5.44) 

o„  -  /)// i,a„.  (.'.45) 

The  appropriate  function-,  for  <!•  as  given  in  table  I.  are  applicable  to  also,  since  the 
symmetry  of  the  applied  loading  is  an  integral  factor  of  the  overall  symmetry.  However,  n  is  not 
necessary  for  the  number  of  parameters  in  a,,  to  match  those  in  a.  The  parameter  series  in  a„ 
is  cut  otf  at  a  point  where  an  adequate  representation  of  the  applied  loading  is  obtained,  bor 
special  applied  load  distributions.  a„  may  contain  only  one  or  a  few  non-zero  parameters.  Kigid 
body  motion  parameters  are  not  applicable  for  an. 


3.7  Stiffness  liquation 

It  has  been  shown  that  the  total  stress  function  for  the  element  can  be  expressed  as: 


-  An  ■  '/'ti- 

(5.46) 

From  this,  it  follows  that 

d  =  d„  do. 

(5.47) 

=  /.,!  a„  /.  not  1 1 . 

(3.4X) 

and 

8  =  Cnaii  ■  titan. 

(5.44) 

so 

an  —  C ii  1  (  ..a.,,  •  C  i,  '8. 

(3.50) 

Also. 

£  -  //uOtu  •  // 1,0,,. 

(5.51) 

so 

a  -  /)£  -  PH,, a„  ■  /)//|,a,,. 

(5.52) 

/)//„«„  /)//„(  t  ,i  'C„au  t  it  '8). 

<  5.55) 

-  />(//„  //„(  it  '(  „)a„  /)//i,Ci,  >8. 

(3.54) 

Now  impose  a  virtual  displacement.  8*. 

d*  -  /. n t "u  '8*. 

(3.55) 

£*  =  //tit'll  '8*. 

(3.56) 

Consider  the  energy  changes  due  to  8*.  Let  q  be  the  vector  of  element  nodal  forces  and  p  the 
vector  of  applied  body  forces. 

(Ill  sternal  work  done  on  the  element: 


by  q.  It  1 4  =-■  (8*)'q.  ( 3.57 ) 

hyp.  It  _  |  (d*)'p</v .  1 3.5X1 

(2)  internal  work : 

Ilf  -  |  (£*)'a(/v.  (3.5d| 

l  or  equilibrium, 

II  :  11,3  =-  II  ,.  ( 5.60) 
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(S*)7^  +  J  (d*)Tpdv  =  J  (e*)roJv,  (3.61) 

(8*)rq  -  J  (HhCh  'S*)T(DH0a0  DHnCu  ‘Cooto-  DHhCh  >6)<A- 

-j*  (in Cm  ‘5*)V-  (3.62) 

(6»)Jq  =  (S*)r(CH~l)rJ  (Hht DHuOLo  —  HliT DH\iCn  -Coat,  -  H,iTDHnCH  '6)‘l\— 

-(8*)T(CH  '*)r  f  LnTpih.  (3.63) 

Since  8*  is  arbitrary, 

q  =  (CH  ')rJ  Hht DH\kI\oq - (C,i  ')tJ  HvJDHndvCn  ‘G>a<>  • 

^(GH  *)rJ  HhtDHh^Ch  '8  (Ch  ■')7'J4  /-HrP</v.  (3.64) 

The  first  two  terms  on  the  right-hand  side  of  (3.64)  are  vectors  of  nodal  forces  related  to  the  initial 
stresses  of  Equation  (3.54).  When  these  terms  were  numerically  evaluated  in  Reference  1  they 
proved  to  be  self  cancelling.  That  result  was  not  able  to  be  proved  analytically  in  Reference  I. 
However,  those  two  terms  are  neglected  in  this  paper  as  having  no  effect.  The  third  term  on  the 
right-hand  side  of  (3.64)  contains  the  stiffness  matrix  of  the  element  (see  (2.2 1 )).  The  last  term  is 
the  standard  formulation  of  the  vector  of  nodal  forces  qi>,  equivalent  to  the  applied  loading 
(see  Ref.  2).  Hence,  (3.64)  reduces  to: 

q=A8-q«.  (3.65) 

Equation  (3.65)  can  be  converted  to  global  cartesian  co-ordinates  using  a  transformation.  T 
(see  (2.22)  to  (2.24)). 

Q  =  *A-Qo.  (3.66) 

Stress  recovery  is  achieved  through  Equation  (3.54)  which  may  be  written  as: 

o=o H  •  v8.  (3.67) 

=  On  t  ,5'A.  (3.68) 

Equation  (3.68)  can  be  evaluated  at  a  number  of  prescribed  points  and  assembled  as  shown 
below. 

f:l  K)  ft-l 


or  =  Oot  i  5 lA  (say). 


3.8  Stiffness  Matrix  Evaluation 

The  main  step  in  evaluating  the  stiffness  matrix  of  an  clement  is  to  evaluate  the  volume 
integral  of  Equation  (2.21).  That  integral  is  given  the  notation: 


imiUh. 


II 


1  ,  r  the  loaded  hole  element,  the  volume  of  integration  is  an  annulus  of  uniform  thickness, 
with  a  circular  inner  boundary  and  a  possibly  irregular,  polygonal  outer  boundary.  In  view  of 
the  complexity  of  the  terms  in  the  //  matrix  (see  Ref.  I ).  it  is  not  feasible  to  evaluate  this  integral 
analytically.  The  method  of  integration  chosen  is  to  approximate  the  outer  polygonal  boundary 
by  a  circular  one.  concentric  with  the  hole,  without  changing  the  volume  of  the  element.  1  his 
allows  the  integration  to  be  performed  between  constant  limits  for  r  and  u.  The  equivalent 
outer  radius  is  given  by: 


(  1  V 

r no  =  'ij  2.  sin(y,.i  y,) 


('.72) 


where  / '  =  0  25  and  0-5  for  the  doubly  and  singly  symmetric  cases  respectively,  and: 

!  I  /\J 

rr:Q  = 


I  /'  1 

wUsi 


n(yi  i  •/,  t  •  sin(yi  y.v ) 


(3.73) 


where /  =  I  0  for  the  unsymmetric  ease:  y,  is  the  angular  co-ordinate  of  node  i.  Now  (3.71 ) 
becomes: 


6 


HTDIIJt>. 


(3.74) 


where  h  is  the  thickness  of  the  clement.  It  was  shown  in  Reference  1  that  the  integial  over  0  is 
readily  performed  analytically,  resulting  in  a  diagonal  type  matrix.  The  integration  over  r  is 
more  complex  and  is  performed  numerically. 


4.  COMPUTER  PROGRAM  DESCRIPTION 

A  computer  program  has  been  written  to  generate  (he  matrices  ami  vectors  A.  Sr.  Q»  and 
am  of  Equations  (3.66)  and  (3.70).  for  a  specified  loaded  hole  element.  The  program  was  written 
in  Fortran  for  use  on  the  PDP-10  computer  of  the  Aeronautical  Research  Laboratories,  in  con¬ 
junction  with  the  general  purpose  finite  clement  analysis  program.  DISMAL.11  The  program, 
entitled  "HOLE",  has  been  entered  in  the  ARL  Computer  Program  Register  and  a  listing  and 
other  details  of  the  program  are  contained  therein. 

The  program  requires  one  input  file  and  generates  three  output  tiles.  The  input  file.  HIN. 
contains  a  geometric  description  of  the  clement,  material  properties,  prescribed  points  for  stress 
evaluation,  and  applied  loading  definition.  For  precise  details  of  the  input  lile  structure  and  a 
listing  of  a  typical  input  file,  see  Appendix  B.  Two  of  the  output  files  are  for  direct  use  by 
DISMAL.  DATA. EXT  is  an  ASCII  file  containing  bookkeeping  information,  and  SPIT  is  a 
binary  file  containing  the  stiffness  and  stress  matrices.  K  and  Sr.  These  files  are  written  in 
APPEND  mode  so  that,  by  re-running  the  program,  the  properties  of  several  loaded  hole  ele¬ 
ments  can  be  strung  together  for  access  by  DISMAL.  The  third  output  lile.  LOADS,  contains 
the  equivalent  nodal  forces  Q»,  the  initial  stresses  o„i.  C  matrix  inversion  parameters,  and  the 
global  stiffness  matrix  A  .  all  in  ASCII  form  and  self  explanatory  format.  Operating  instructions 
for  the  loaded  hole  element  program  in  conjunction  with  DISMAL  are  given  in  Appendix  C. 

The  program  caters  for  three  modes  of  specifying  the  applied  loading  (see  Appendix  A.) 
Firstly.  Fourier  coefficients  can  be  input  directly.  The  program  user  must  c  sure  that  the  input 
coefficients  conform  to  the  allowable  functions  appropriate  to  the  particular  loaded  hole  element 
variant  (see  Table  I).  Secondly,  point  loads  at  various  angular  co-ordinates  can  be  specified. 
The  program  automatically  allows  for  symmetry  in  calculating  the  Fourier  coefficients  for 
point  loads,  by  excluding  inappropriate  coefficients.  Finally,  for  general  distributed  loading,  the 
program  reads  ordinate  values  at  an  odd  number  of  equally  distributed  points  aiound  the  hole 
boundary.  Again,  the  program  automatically  enforces  symmetry.  For  the  double  and  singly 
sy  mmetric  variants  a  reflected  set  of  ordinates  is  generated  from  the  input  ordinates.  The  program 
obtains  the  Discrete  Fourier  Transform  of  the  ordinate  values  by  calling  a  subroutine  entitled 
C06AAF  from  the  NAG  library  of  subroutines  available  on  the  PDP-10  computer  at  ARL. 
For  the  doubly  symmetric  variant  the  frequencies  of  the  computed  Fourier  coefficients  are 
doubled,  to  allow  for  the  fact  that  the  input  and  reflected  ordinate  values  cover  a  range  *  rather 
than  2n.  f  or  representing  the  applied  loads,  it  has  proved  necessary  to  use  «  -  8  as  a  limiting 
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frequency,  as  inclusion  of  higher  frequency  coefficients  introduces  large  numerical  errors  into 
the  calculation  of  the  initial  stress  ou.  The  program  caters  for  the  specification  of  no  applied 
loading,  in  which  case  it  functions  as  an  unloaded  hole  element.  No  equivalent  nodal  forces  or 
initial  stresses  are  then  output. 

The  integration  over  r  of  Equation  (3.74)  is  performed  numerically,  using  Simpson's  rule 
with  100  points.  The  inversion  of  the  matrix  Cu  of  Equations  (3.49)  and  (3.64)  is  performed  by 
a  subroutine.  GENINV.7 

The  points  at  which  the  vector  oo  and  the  matrix  S  of  (3.67)  are  evaluated,  to  construct  the 
vector  Oo  r  and  the  matrix  St  of  (3.70).  consist  of  radial  projections  of  the  nodal  points  onto  the 
hole  boundary,  other  prescribed  points  on  the  hole  boundary,  and  general  prescribed  points. 
For  points  on  (he  hole  boundary,  only  the  tangential  stress  component  is  included  in  Out  and  St. 
as  the  radial  and  shear  stress  components  are  known  from  the  applied  loading.  Co-ordinates 
of  prescribed  points  are  defined  in  the  input  file  to  the  program. 

In  principle,  it  is  not  necessary  to  have  a  uniform  distribution  of  nodes  around  the  circum¬ 
scribing  circle  of  the  element.  However,  unexplained  numerical  difficulties  have  been 
encountered  when  using  the  program  with  non-uniform  distributions  of  nodes  (see  Ref.  I). 
Therefore,  the  use  of  the  current  loaded  hole  clement  program  is  restricted  to  uniform  distribu¬ 
tions  only.  The  aspect  ratio  of  the  element,  n  m.  is  unrestricted  except  to  be  greater  than  unity. 
However,  results  of  example  analyses  indicate  that  this  ratio  should  be  in  the  range  I  -5  to  3  0 
for  best  accuracy  (see  Ref.  I).  Other  restrictions  on  the  input  data,  necessary  for  satisfactory 
operation  of  the  program,  are  listed  in  Appendix  B. 

5.  EXAMPLE  ANALYSES 
S.l  Problems  and  Results 

Many  example  analyses  xverc  described  in  Reference  I.  Some  of  those  are  included  in  this 
section  in  order  to  demonstrate  the  use  of  the  element  and  the  accuracy  to  be  expected  from  it. 
The  first  example  is  a  rectangular  plate  with  a  central  circular  hole,  subject  to  uniform  longi¬ 
tudinal  tension  (see  Fig.  3).  This  is  a  doubly  symmetric  problem  and  the  mesh  used  is  shown  in 
Figure  3.  A  doubly  symmetric  hole  element  with  eight  uniformly  spaced  nodes  was  used.  The 
maximum  stress  concentration  around  the  hole  boundary  is  of  interest,  i.e.  the  tangential  stress 
n„„  at  r  =  ro.  0=0  divided  by  the  uniform  applied  stress.  Accurate  values  for  this  stress  con¬ 
centration  factor,  accounting  for  the  finite  width  of  the  plate  but  not  the  finite  length,  are 
available  in  Reference  X.  This  problem  was  analysed  for  a  range  of  values  of  r, i,  keeping  the 
element  circumscribing  radius,  constant.  Zero  applied  loading  was  specified  for  the  hole 
element.  The  results  are  shown  in  Table  3. 

TABLE  3 


Stress  C  oncentration  for  Hole  in  a  Rectangular  Plate 


n  i 

n  /  ii 

(f,w/)ina.\  fT 

ts.c.r.  )„-,■* 

error 

1  0 

4  0 

3-56 

3  02 

17-8 

1  -25 

3  2 

3  02 

3-04 

0-7 

1-5 

2-67 

3-05 

3  05 

0  0 

2  0 

2  0 

3  12 

3  10 

0  7 

2-5 

1  -6 

3-20 

3-16 

1-2 

3  0 

1-33 

3  30 

3  -  24 

1  9 

3  5 

1  14 

3-44 

3  34 

3  0 

*  Accurate  value  of  stress  concentration  factor.  "■  obtained 

from  Reference  X. 


The  performance  of  the  singly  symmetric  variant  of  the  loaded  hole  element  is  demonstrated 
by  the  problem  of  an  oll'sct  circular  hole  in  a  circular  disc,  subject  to  uniform  internal  pressure. 
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FIG  3:  HOLE  IN  RECTANGULAR  PLATE 


The  mesh  used  is  shown  in  Figure  4.  I'niform  internal  pressure  of  I  u  ■  It)'  is  •■pecilie.1  h. 
inputting  the  frequency  coeflieient  aa „  =  2  0  <I0’  with  n  -  ()  (see  \p|vn.iiv  Hi  V.  urate 

values  for  the  stress  concentration  factor  at  points  A  and  B  of  Figure  4  lute  Iven  ohla.ucu  front 
Reference  8.  Values  derited  from  the  analysis  of  the  ntesli  of  Figure  4  hate  been  soniparcd  to 
them  in  Table  4. 

TABLE  4 

Stress  Concentration  for  an  Offset  Hole 

/>»  (s.c.f.)aiv  I  error 

| 

2-11  217  I  24 

1  68  1  65  |  14  ' 

|  I 

An  example  problem  with  a  broad  spectrum  of  frequency  components  in  the  loading 
applied  to  the  hole,  is  a  riveted  joint  with  load  transfer  through  the  rivets.  The  load  distribution 
is  of  cosine  form  over  the  area  of  contact  between  the  rivet  and  the  little.  I  his  problem  is  detailed 
in  Figure  5.  For  the  riveted  joint  problem  the  applied  loading  was  input  as  ordinate  values  and 
the  frequency  components  established  by  taking  the  1)FT.  The  tangential  stress  was  calculated 
at  many  points  around  the  hole  boundary  and  the  maximum  stress  concentration  factor. 

was  found  lo  be  2- 15  at  0  =  58  .  This  is  5 •  7"„  in  error  with  an  accurate  value  ol 

2 ■  28  at  0  =  60  given  in  Reference  8. 

The  versatility  of  the  loaded  hole  element  is  demonstrated  by  the  analysis  of  two  holes,  of 
generally  unequal  diameter,  in  a  large  plate  (clleciively  infinite).  Ihe  ineslt  used  ts  shown  in 
Figure  6.  For  each  run  of  the  problem  two  load  eases  were  analysed,  v  i/.  each  hole  u,  turn  subject 
to  uniform  internal  pressure  />»  with  the  other  hole  unloaded.  From  urn  to  run.  the  radii  of  t lie- 
holes  were  varied  without  altering  the  mesh  and  keeping  n  r„  in  the  langc  I  5  to  t  n.  I  or  either 
load  ease  of  any  run.  the  maximum  tangential  stress  occurs  at  a  point  C  on  the  loaded  hole 
surface,  displaced  by  an  angle  (V  from  the  line  joining  the  hole  centres.  I  he  'tress  concentrations 
at  point  A  on  the  loaded  hole  boundary  and  point  B  on  the  unloaded  hole  boundary  ate  also  of 
interest.  Points  A.  B  and  C'  are  illustrated  in  Figure  7.  |  he  results  front  all  runs  are  given  in 
Table  5  where  the  subscript  I  refers  to  the  loaded  hole  and  the  subscript  2  lo  the  unloaded  hole. 
Ihe  slress  concentrations  at  point  C  in  fable  5  have  been  plotted  in  I  igurc  7.  I  or  the  combined 
loading  ease  of  both  holes  under  equal  pressure  when  rm  -  inj.  the  stress  concentration  at  points 
A  and  B  can  be  found  by  adding  the  values  in  columns  A  and  B  of  ruble  5.  I  hew  combined 
values  are  compared  with  accurate  values  from  Reference  8  in  Table  6. 

5.2  Discussion  of  Results 

Ihe  example  analyses  outlined  in  the  preceding  section  demonstrate  the  accuracy  and 
utility  of  the  loaded  hole  element.  The  lirst  example,  a  central  hole  in  a  rectangular  plate,  indiv  atc-s 
that,  for  best  accuracy,  the  aspect  ratio  of  the  element,  n  r».  should  he  kept  in  the  range  I  -5 
to  2  0.  For  that  problem,  errors  were  then  less  than  2",  .  I  he  analysis  of  the  riveted  joint  shows 
that  the  clement  can  produced  accurate  stress  concentration  values,  even  for  complex  applied 
load  distributions.  The  mesh  for  that  problem,  shown  in  figure  5.  uses  several  six-nodevl  tri¬ 
angular  elements,  two  of  which  adjoin  the  loaded  hole  element.  Ihe  two  mid-side  nodes  along 
the  junction  arc  not  nodes  of  the  hole  clement.  I  nconnectcd  nodes  are  normally  undesirable 
as  introducing  inter-element  displacement  incompatibility.  Howexer.  xvith  standard  element 
connections  to  Ihe  loaded  hole  element,  displacement  incompatibility  is  inevitable  anyway, 
and  using  unconnected  nodes  will  not  necessarily  worsen  it.  Ihe  analvsis  ol  two  adjacent  holes 
in  an  infinite  plate  shows  the  flexibility  of  the  loaded  hole  element.  I  sing  only  one  mcslt.  broad 
ranges  of  the  basic  parameters  ol  the  problem  were  spanned,  with  errors  less  than  3",..  An 
alternate  solution  of  this  problem  is  unknown,  and  the  stress  concentration  values  in  lablc  5 
arc  then  useful  basic  data. 
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FIG  6:  TWO  HOLES  IN  AN  INFINITE  PLATE 


TABLE  5 


Stress  Concentration  for  Two  Holes  in  an  Infinite  Plate 


rni 

ny> 

r  hi  /r 02 

1 

2/oi  'h 

j  <?„„  pO 

i 

I 

"<■ 

(deg) 

A 

B 

c 

10 

1  0 

10 

0-333 

1  -007 

0-  163 

1  -014 

70 

1  -0 

1-5 

0-667 

0-333 

1  -000 

0-199 

1  027 

60 

15 

10 

1-5 

0-5 

0  979 

0-374 

1-031 

' 

80 

1  0 

2  0 

0-5 

0  333 

0-989 

0  249 

1  -050 

60 

20 

10 

2-0 

0-667 

0  963 

0  664 

1  -042 

80 

1-5 

1-5 

1-0 

0-5 

0  991 

0-460 

1  062 

60 

1-5 

20 

0-75 

0-5 

0  953 

0-589 

11 18 

60 

2  0 

1  -5 

1  •  333 

0-667 

0  906 

0-835 

1  •  103 

40 

O 

ri 

20 

1  -0 

0-667 

0-866 

1  093 

I  -254 

40 

TABLE  6 


Stress  Concentration  for  Two  Equal  Holes 
in  an  Infinite  Plate 


'ul  =  r  n_’ 

h  2/-,h 

{ ^ lift  /'ll)p«tini  A 

error 

Hole  element 

Reference  8 

1  -o 

3 

1-170 

1-153 

1-5 

I 

1-5 

2 

1  -451 

1  -410 

i 

2-9 

2-0 

1-5 

1  -959 

1  900 

3-1 

The  errors  listed  in  Table  ft  can  be  considered  as  upper  limits  on  the  possible  errors  in 
Table  5. 

6.  DISCISSION 

The  loaded  hole  element  has  many  advantages  o\et  an  equivalent  mesh  of  standard  elements, 
firstly,  stresses  can  be  obtained  accurately  at  any  point  within  the  element.  Specifically,  the 
process  of  matching  boundary  stresses  with  applied  loads  ensures  accurate  stresses,  at  or  near 
the  boundary.  Standard  elements  cannot  produce  boundary  stresses,  no  matter  how  much  the 
mesh  is  relined.  In  order  to  obtain  boundary  stresses  from  them,  it  is  necessary  to  extrapolate 
to  the  boundary,  and  in  the  case  of  a  hole  boundary,  the  extrapolation  is  into  a  region  of  high 
stress  gradients.  When  using  the  loaded  hole  element,  both  preparation  and  computing  time  are 
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cm  significant!)  from  an  equivalent  anal) sis  with  a  standard  element  mesh.  This  is  due  to  a  reduc¬ 
tion  in  the  number  of  elements  and  nodes.  The  displacement  incompatibility  of  the  loaded  hole 
element  with  adjoining  standard  elements  did  not  prove  to  he  a  significant  factor.  Bounds  on  the 
accuracy  of  the  solutions  to  example  problems  were  not  established,  due  to  the  oscillatory 
nature  of  the  convergence.  However,  for  most  engineering  applications,  all  that  is  required  is 
a  knowledge  of  the  "v  peeled  accuracy  of  the  analysis,  and  experience  in  that  regard  can  be  gained 
from  comparisons  of  analyses  using  the  loaded  hole  element  with  theoretical  solutions  (see 
Section  5). 

As  mentioned  in  the  introduction  to  this  report,  the  hole  element  was  developed  as  a  stage 
in  the  development  of  an  element  capable  of  representing  a  warhead  fragment  damage  site  in  an 
aircraft  structure.  Such  a  damage  site  can  be  idealised  as  a  hole  with  cracks  radiating  from  it. 
As  vet.  endeavours  to  extend  the  hybrid  element  method  to  model  a  hole  with  one  or  two  long 
cracks  radialin  fn  m  it  have  met  with  numerical  difficulties,  and  this  goal  has  had  to  be  left  for 
later  study. 

7.  OOM'l.l  SION 

The  satisfactory  operation  of  the  three  variants  of  the  loaded  hole  element  has  been  con¬ 
ed  for  uniform  distributions  of  4  to  10  nodes.  Adequate  operating  instructions  for  the  loaded 
h.  le  element  program  are  contained  in  Appendices  B  and  C'  to  this  report.  A  listing  of  the  pro¬ 
gram  is  contained  in  the  ARL  Computer  Program  Register. 
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APPENDIX  A 

Determination  of  Fourier  Coefficients 


The  applied  loads  are  represented  as: 


Vs 

a„(l>)  =  CT„(,  2  -  X  l"an  COS  nil  7  CTan  sin  III1 ; 

« -  t 

's 

T»(I0  =  ra(,  /2  •  V  |ra„  sin  /ffi  -  Ta„  COS  /id; 

n  ~  i 


Where 


(A.I) 
(A. 2) 


1  C'-” 

"an  =  1  "alffi  COS  llOl/l) 

J  0 

(A. }) 

|  fin 

"an  =  |  "affi)  sin  IlflllO 

"J  0 

(A. 4) 

ra„  —  f  Ta(d)  sin  nfldO 
irj  0 

1  /*-" 

(A. 5) 

ra„  =  ra(0)  cos  IllklO 

"  J  » 

(A. 6) 

.  ra„.  ra„  can  be  determined 

in  three  ways  depending  on  the 

method  of  specifying  the  applied  loading.  oa(<t). 


(1)  If  the  applied  stress  distributions,  ua(ffi  and  ra(/>)  are  trigonometric  as  sin  nO  or  cos  nfl. 
then  the  coefficients  are  known  directly.  Also  if  the  distributions  oa('0  and  ra(U)  are 
known  functions,  such  that  the  integrals  of  equations  (A. 3)  to  (A. 6)  can  be  easily 
evaluated  analytically,  then  the  coefficients  can  be  evaluated  external  to  the  program 
and  input  to  it. 

(2)  As  a  special  case  of  known  functions,  if  point  loads  are  applied  then  od{H)  and  ?d(0) 
are  delta  functions  of  the  form: 


oa(0)  =  0f) 

fob 

(A. 7) 

S 

~  Of) 

r^b 

(A. 8) 

Where  Pd i  and  Sm  arc  radial  and  shear  loads  applied  at  0<.  ru  is  the  hole  radius,  and  h  the  element 
thickness.  Then  Equations  (A. 3)  to  (A. 6)  become: 


P at 

cos  nth 

niob 

(A.9) 

P at  . 

.  sin  nth 

TTfub 

(A.  10) 

Sm 

.  sin  nth 
nruo 

(A.I  1) 

*Sa  t 

.  cos  n*h 

nfub 

(A. 12) 

I 

t 
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(3)  If  n^O)  and  ra(d)  are  complex  or  unknown  functions,  they  can  be  input  as  ordinate 
\alues  at  A’„  discrete  values  of  0.  The  points  are  equally  spaced  and  labelled  0  •  X,,  I 
where  .V,,  is  an  even  number.  The  distributions  0^(0)  and  ra(/q  are  then  represented  as 
Fourier  series  by  taking  the  Discrete  Fourier  Transform  of  these  points.  The  DF'T  is 
given  by: 


(Tan  i  *  <7o  4  (  ;V 

<T*{0)  =  *  £  {°an  COS  tlQ  *-<Tan  sin  tlO\  *-  COS-  ”0 

2  n  =  l  2(2 

,“i  ! 

T»(h)  =  T“°  r  Z  { Ta«  cos/ifl-Ta„  sinnUJ  f  Tat.'’  cosj  S' 

2(2 

! 

(A. 14) 

Where 

; 

2  ■'>  -  i 

°ai>  =  Z  COS(2»imr//Vp) 

'*P  »i  =  ' 

1 

(A. 15) 

2  »y  i 

=  ,  Z  sin(2/WW/;V  ) 

Apm  -  | 

(A. 16)  | 

2  v„  -  i 

^  Z  Ta(Wm)  COS(2lHHir/iV  ) 

A  =  o 

(A. 17) 

2  N>  =  1 

Ta*i  =  v  Z  Ta(^/«)  sin(2m/;WAp) 

m  =  1 

(A.  IX)  | 

Equations  (A. 15)  to  (A. 18)  are  numerical  evaluations  of  equations  (A.3)  to  (A. 6).  The 
coefficients  aa„  .  cra„  .  ra„  ,  ra„  are  evaluated  using  a  DFT  package  program.  The 

DFT  can  be  put  in  the  form  of  Equations  (A.I)  and  (A. 2)  by  using  the  following 

relationships. 


°an 

= 

(A. 19) 

a*  n 

= 

(A.20) 

Ta« 

= 

r  an 

(A. 21 ) 

Tan 

= 

~  T&ll 

(A. 22) 

As 

= 

'V2 

(A. 23) 

= 

0-5<ra|.v, 

(A. 24) 

Ta.Vs 

= 

0-5rai.v, 

(A. 25) 
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APPENDIX  B 


Input  Data  File 

The  input  data  file  is  entitled  HIN  and  contains  all  parameters  necessary  to  define  the 
clement  geometry  and  material  properties,  the  stress  evaluation  points,  and  the  applied  loading. 
The  layout  of  the  parameters  in  file  HIN  is  given  below. 

///.V 

A>.  r  i.  b.  f 
V 

}'<l).  y<2) . y(\  ) 

. 

•Vo.  I  n.  .Vi.  V, 

A  It 

l‘i i(  I ).  hut 2) . t'n(.V),)  lOmit  if  .V,,  =  O; 

\,i 

/«;( I ).  ",  A  1 1.  r,A  2).  <»,;(2) . lOmit  if  .V(i  =  0; 

P  1=  TRIG  | 

At 

II.  n aw.  (J an  .  r.t((  ,  t;,m 

ii.  (TU„.  ■  ■  ■  fete,  for  A't  lines ; 

P  I  =  POINT  | 

■V,..  As 

"i-  Pat.  Sat 

l'i.  P*i.  .  .  .  {etc.  for  N r  lines; 

P  {  =  NUMER | 

Mi 

rta(U.  Oa(2).  oa(3) . <ra(A'l-d) 

ra(  1 1.  Ta(2l.  Ta(3) . ra(tVpn) 

P  {  =  END  I 

The  statements  enclosed  in  curly  brackets  are  comments  and  not  part  of  the  file.  The  file  segments 
delineated  by  square  brackets  can  be  placed  in  any  order  and  any  segment  type  can  be  repeated 
or  omitted.  For  no  applied  loading  all  three  segments  should  be  omitted.  All  quantities  in  the 
tile  are  in  list  directed  formal.  Where  practicable,  the  parameter  symbols  in  the  above  file  layout 
correspond  with  symbols  used  in  the  main  body  of  this  report,  rather  than  with  the  variable- 
names  used  in  the  program.  For  clarity,  all  symbols  are  defined  in  Table  Bl  and  their  limitations 
prescribed.  All  parameters  referred  to  in  Table  Bl  as  "numbers"  are  integers.  Other  unspecified 
parameters  are  real  numbers.  All  angles  are  input  in  degrees.  Other  quantities  are  input  in  any 
self-consistent  set  of  units.  The  program  was  developed  for  anticlockwise  node  numbering  and 
its  successful  operation  for  clockwise  numbering  has  not  been  confirmed.  For  applied  loading, 
the  sign  convention  is  positive  for  inward  radial  loads  and  anticlockwise  shear  loads.  Naturally, 
all  limits  specified  in  Table  Bl  ire  in  addition  to  inherent  sensibility  limits,  e.g.  it  would  be  non¬ 
sense  to  specify  a  negative  element  thickness,  h. 


) 


As  an  example,  the  data  file  HIN  for  a  doubly  symmetric  loaded  hole  element  with  10  nodes 
is  listed  below.  The  hole  radius  is  2  0m.  the  element  circumscribing  radius  is  4  0  m  and  the 
element  thickness  0  002  m.  Stresses  are  to  be  evaluated  at  three  hole  boundary  points  only. 
Applied  loading  consists  of  a  point  load  superimposed  on  a  distributed  load. 

HIN 

20.  4  0.  0  002,  0-25 
10 

0  0,  10  0,  20  0,  30  0.  40  0.  50  0,  60  0.  70  0,  80  0,  90  0 
7  I0E07.  0-3 
0  0,  0  0.  I  0.  0  0 

3 

5  0,  15  0,  25  0 
0 

POINT 

I.  8 

45  0,  1000  0.  0  0 
Nl'MER 

4 

I  0E05,  I  2E05.  I  4E05.  I -3E05,  I1E05 
0  0.  I  -5E05.  I  -6L05,  l  5E05,  0  0 
END 


TABLE  B1 

Input  Parameter  Definitions 


1 

Symbol 

Definition 

j 

1 

Limits 

n> 

Hole  radius 

In  >  0 

n 

Element  circumscribing  radius 

>i  >  rn  and  I  S  <  ii  <  3  0  for 
best  accuracy 

b 

Element  thickness 

Validity  of  plane  stress  assumption 

/ 

Element  fraction  of  a  full  circle 

=  0-25  for  doubly  symmetric 
=  0-5  for  singly  symmetric 
=  10  for  unsymmetric 

.V 

Number  of  nodes 

3  ^  AA  10 

yU) 

Angular  co-ordinate  of  node  i 

■/(it  >  all  y(j).  for  j  <  i 

7(A)-  7(1)  =  360 f.  for  /*  10 
<  360 f  for /  =  10 

E 

Young's  modulus 

None 

V 

Poisson's  ratio 

None 

•V».  )  « 

Globul  cartesian  co-  ordinate  of  the  hole 

centre 

None 

•V,.  )i 

Global  cartesian  co-ordinate  of  node  1 

None 

An 

Number  of  hole  boundary  points  for 
stress  evaluation 

/V  •  A’u  -  3/Yt;  ^  30 

fhtU) 

Angular  co-ordinate  of  boundary  stress 
evaluation  point  i 

None 

An 

Number  of  general  points  for  stress 
evaluation 

A  •  A„  -f-  3 A’u  $  30 

rr.(i).  Oc.U) 

Local  polar  co-ordinates  of  general 
stress  evaluation  point  i 

None 

P 

A  flag  to  specify  the  type  of  loading 

TRIG.  POINT.  NUMER. 

END 

At 

Number  of  frequencies  for  which 
Fourier  coefficients  are  input 

None 

n 

Frequency  parameter 

Integer.  0  s*  n  Si  16 

^itw-  ^ari 

Coefficients  of  cos  nO  and  sin  nil  for 
applied  normal  stress 

These  functions  are  referenced  to  node 

1  as  n  -  0 

ran  .  tbh 

Coefficients  of  cos  uf  and  sin  nO  for 
applied  shear  stress 

These  functions  are  referenced  to  node 

1  as  0  =  0 

An 

Number  of  positions  at  which  point 
loads  are  applied 

None 

As 

Cut-off  frequency  for  Fourier  represen¬ 
tation  of  point  loads 

Integer.  0  $  A's  <  16 

Ih 

Angular  co-ordinate  of  point  load 
position  i 

None 

Pm.  Sm 

Radial  and  shear  point  load  compon¬ 
ents  at  position  i 

None 

Mi 

Parameter  defining  the  number  of  input 
ordinate  values.  A’fo 

Even  integer  <  8.  10.  12  for /  =0-25, 
0-5.  10 

A  FO 

Number  of  input  ordinate  values 

Af„  =  2-"i  -  t  1 

"a(0,  Ta(l) 

The  /th  normal  and  shear  stress  ordinate 
values 

Ordinate  values  are  for  equally  spaced 
points  around  hole  boundary.  For 
f  —  0-25.  0-5.  points  1  and  Afii  lie 
at  nodes  1  and  A  respectively,  and 
rati)  =  t»(A'fh)  =0.  For  f  —  10. 
points  1  and  A’fo  lie  at  node  1  and 
u«(l)  =  <)„(  A'fu)  and  T»(l)  =  T,(A'fo) 

APPENDIX  C 


Operating  Instructions 

The  step-by-step  operations  for  using  the  loaded  hole  element  program.  HOLE,  in  con¬ 
junction  with  the  general  finite  element  analysis  suite  of  programs.  DISMAL,  are  listed  below. 

(1)  Draw  up  the  mesh  for  the  structure,  including  any  loaded  hole  elements,  and  define  the 
applied  loads. 

(2)  Prepare  the  data  file.  DATA,  for  DISMAL  (see  Ref.  6).  excluding  the  applied  loading 
data  at  this  stage.  Loaded  hole  elements  are  specified  in  DATA  as  special  elements 
(see  Ref.  6). 

(?)  For  each  loaded  hole  element  (in  the  order  of  specification  in  DATA)  prepare  a  data 
file,  HIN  (see  Appendix  B).  run  the  program  HOLE,  and  obtain  a  listing  of  CHECK.  LST. 

(4)  Complete  the  applied  loading  section  in  DATA  including  the  equivalent  nodal  forces 
obtained  from  the  listings  of  CHECK. LST. 

(5)  Run  the  DISMAL  programs  (see  Ref.  6). 

(6)  Obtain  the  stresses  for  the  loaded  hole  elements  from  the  DISMAL  output  file 
PRINT3.LST.  These  stresses  will  be  output  in  the  order  of:  tangential  stress  for  nodal 
point  projections  followed  by  specified  boundary  stress  evaluation  points;  then  followed 
by  sets  of  radial,  tangential  and  shear  stress  for  general  stress  evaluation  points.  Add 
these  stresses  to  the  corresponding  initial  stresses  obtained  from  the  listings  of 
CHECK. LST. 

The  above  instructions  apply  to  the  case  w  here  only  one  load  case  is  specified  in  DATA.  If  more 
than  one  load  case  is  to  be  specified,  involving  different  loading  on  the  holes,  it  will  be  necessary 
to  establish  equivalent  nodal  forces  and  initial  stresses  by  preliminary  runs  of  HOLE,  and  then 
establish  the  correctly  structured  files  DATA. EXT  and  SPIT  by  an  appropriate  sequence  of 
runs  of  HOLE.  The  loaded  hole  element  can  be  used  in  conjunction  with  other  special  elements, 
providing  the  correct  appending  of  files  DATA. EXT  and  SPIT  is  monitored. 
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